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Diffuse UV light associated with the Spiderweb Galaxy: evidence for 
in-situ star formation outside galaxies 
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ABSTRACT 

We present detailed images of diffuse UV intergalactic light (IGL), situated in a 60 kpc halo 
that surrounds the radio galaxy MRC 1 138-262 at z=2. We discuss the nature of the IGL and 
rule out faint cluster galaxies, nebular continuum emission, synchrotron, inverse Compton, 
synchrotron self-Compton emission and scattering of galactic stellar light as possible sources 
of the IGL. Dust scattered quasar light is an unlikely possibility that cannot be ruled out en- 
tirely. We conclude that the source of the IGL is most likely to be a young stellar population 
distributed in a halo encompassing the radio and satellite galaxies, undergoing star formation 
at a rate greater than 57±8M @ yr . Within 70 kpc of the radio core, approximately 44% 
of the star formation that is traced by UV light occurs in this diffuse mode. The average UV 
colour of the IGL is bluer than the average galaxy colour, and there is a trend for the IGL to be- 
come bluer with increasing radius from the radio galaxy. Both the galaxies and the IGL show 
a UV colour-surface brightness relation which can be obtained by variations in either stellar 
population age or extinction. These observations show a different, but potentially important 
mode of star formation, that is diffuse in nature. Star formation, as traced by UV light, occurs 
in two modes in the high redshift universe: one in the usual Lyman break galaxy clump-like 
mode on kpc scales, and the other in a diffuse mode over a large region surrounding massive 
growing galaxies. Such a mode of star formation can easily be missed by high angular resolu- 
tion observations that are well suited for detecting high surface brightness compact galaxies. 
Extrapolating from these results, it is possible that a significant amount of star formation oc- 
curs in large extended regions within the halos of the most massive galaxies forming at high 
redshift. 

Key words: galaxies:individual:MRC 1 138-262 - galaxies:haloes - galaxies:high-redshift - 
galaxies: elliptical lenticular, cD - galaxies:clusters:general. 



1 INTRODUCTION 

Distant powerful radio galaxies are important probes of the for- 
mation and evolution of massive galaxies, because they are among 
the most luminous and massive galaxies known in the early Uni- 
ver se, and are likely progenitors of dominant cluster galaxies 
(e.glPentericci et alJl997MVillar-Martm et al.l200rj : lsevmour et all 
. They are generally embed ded in giant (cD-sized) ionized 



gas halos (e.glReuland et al.l2003 ) and surrounded by ga laxy over- 
densities (Pent ericci et al.l 12000: Ven emans et al. 2007). With ra- 
dio lifetimes (few x 10 7 yr) being much smaller than cosmological 
timescales, the statistics are consistent with every dominant clus- 
ter galaxy having gone through a luminous radio phase during its 
evolution. 

The deepest image of a distant radio galaxy with the Hub- 
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ble Space Telescope (HST) is a recent Advanced Cam era for Sur- 
veys ( ACS) observation of MRC 1138-262 at z= 2.156 ( Mile y et al.l 

2006) . This system has been dubbed the Spiderweb galaxy as its 
complex morphology resembles a spider's web. The ACS image 
shown in Fig. 1 shows tens of UV bright galaxies surrounding a 
central galaxy (marked with a white cross). As the radio core is 
cospatial with the central galaxy we refer to this central galaxy as 
the radio galaxy, and all the surrounding galaxies are referred to as 
the satellite galaxies. All of these galaxies may merge together to 
form a single massive galaxy at z=0, therefore the whole complex 
system is referred to as the Spiderweb system This system is one of 
the most intensive ly studied distant radio galaxies l lPentericci et al.l 
1997, 1998, 2000) and exhibits the properties expected for a pro- 
genitor of a dominant cluster galaxy. The infrared l uminosity pro- 
vides an upper limit to the stellar mass of 10 12 M© ^Seymour et al.l 

2007) , so it is one of the most massive galaxies known at z>2. The 
host galaxy is surrounded by a giant Lya halo and the high rotation 
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measure of the radio source means that it is embedded in a dens e 
medium with an ordered magnetic field dPentericci et al.l [l997l) . 
The radio galaxy is associated with a >3 Mpc-si zed structure of 
galax ies with an estimated mass > 4 x 1O 14 M dVenemans et al.l 
1 20071) . indicating that it is the predecessor of a local rich cluster. 
Mechanical feedback from the active galactic nucleus (AGN) is 
sufficient to expel significant fra ctions of the interstel lar medium 
of the massive gas-rich galaxy dNesvadba etai] [2006). Therefore 
the AGN may quench the star formation and allow the radio galaxy 
to evolve onto the red seque nce as suggested b y popular models 
of massive galaxy formation dCroton et alj|2006l) . The study of the 
formation of the most massive galaxies through these and other 
mechanisms has so far mainly been limited to models and sim- 
ulatio ns (e.g |Pubinski| [i998; Gao et al. 2004; |Pe Lucia & Blaizotl 
l2007h . Our aim is to use these observations of a high-redshift radio 
galaxy to study brightest cluster galaxy formation. 

In this work we present detailed images of widespread UV 
intergalactic light (IGL) situated in a 60 kpc halo that lies between 
the radio galaxy and the UV bright satellite galaxies. We examine 
the possible origins of this light and present in-situ star formation 
as the most probable. In section [2] we describe the observations 
and data reduction, in section [3] we describe the distribution and 
colour of the IGL, and discuss possible origins in section|4] Section 
[5]discusses the implications of widespread star formation in a halo 
around a massive forming brightest cluster galaxy. Thr oughout this 
work we use H = 71, Q M = 0.27, and Q. A = 0.73 dSpergel et al] 
2003). All magnitudes are AB magnitudes. At a redshift of 2.156, 
the linear scale is 8.4 kpc/". 



2 METHOD 
2.1 Observations 

The imaging data used in this work was obtained with the ACS 
on the HST. The total exposure time in the F475W (#475) fil- 
ter was 9 orbits (20670 s) and 10 orbits (23004 s) in the F814W 
(hu) filter. The data were reduced using the ACS pipelin e sci- 
ence investigation software (Apsis; iBlakeslee et al] 120031) . The 
5cr depth of the 5475 and /gi4 images is 27.5 mag arc sec _ 2 and 
27.0 mag ar csec -2 respective ly. Further discussion of the ACS data 
is given in Miley et al. ( 20o3)- For our analysis we extract a region 
of 15"x 18.5" around the central radio galaxy which avoids any ob- 
vious bright foreground stars. The combined £475 + ^814 image of 
this region is presented in Fig. [7] The total column density of our 
galaxy in the direction of the radio galaxy is 4.6 x 10 20 cm~ 2 which 
results in a galactic reddening of E (B — V) = 0.079. 

The UV spectral data were obtained using the multi object 
spectrograph (MOS) on FORS1 in ESO program 64.O-0 134(D). 
These data were obtained as a by-product of the confirmation of 
the overdensity o f Lya emitters in the surrounding protocluster 
dKurk et al] 2000). The position of the slit on the sky is shown in 
Fig. Q] The spectra were obtained with the 600B grism, and a 1" 
wide slit, giving a dispersion of 1.2 ApixeP 1 and a pixel scale of 
0.2" in the spatial dimension. The resolution of the spectra was ap- 
proximately 5 A /400km s . The total exposure time was 6 hours 
(divided in 6 separate exposures of 1 hour each). The seeing con- 
ditions varied between 0.7" to 1.5". The data were reduced using 
standard IRAF tasks, which included: bias subtraction, flat field 
removal, and background subtraction. A relative flux calibration 
was achieved using the spectrophotometric standard star GD108, 
whilst wavelength calibration was performed using the arc-lamp 



exposures (He and HgCd lamps). The accuracy of the wavelength 
calibration is estimated to be 0.2 A. 



2.2 Emission lines within the g475 and Igu bands 

The £475 and 7^14 wavelength ranges both encompass a num- 
ber of emission lines, notably CIVA1549 and HeIlAl640 in the 
g 475 , and Cll]A2326, [NeIV]A2424,2426, and MgIlA2798 in the 
/814. We do not have a spectrum of the emission-line halo sur- 
rounding MRC 1138-262 that covers the lines mentioned redward 
of HeIlAl640. The Cll], [Neiv] doublet, and Mgll ar e typically 
weak: the composite radio galaxy spectrum presented bv lMcCarthvl 
( 1993) shows that these lines have equivalent widths of 19A, 20A 
and 19A respectively and therefore make up only 7 per cent of the 
flux within /814. We therefore assume that the emission line con- 
tribution to the 7gi4 flux is negligible. McCarthjj d 19930 give the 
combined equivalent width of the Civ and Hell lines of a compos- 
ite radio galaxy to be 136A which results in 34 per cent of the flux 
in the £475. We therefore use spectra of the Lya halo of MRC 1 138- 
262 to estimate the contribution these lines have to the #475 flux. 

To calculate the exact contribution of the emission lines to the 
g475 flux we use the FORS 1 long-slit spectrum to examine the UV 
spectra of 3 distinct regions in the field-of-view: the radio galaxy; a 
nearby satellite galaxy (labeled galaxy A in Fig.[TJ; and the region 
in the remainder of the slit (both East and West of the radio galaxy) 
that did not have any peaked continuum in the spatial direction but 
had low level continuum and Lya emission - labelled the diffuse 
continuum region. The spectra are displayed in Fig.|2]and show that 
the emission line contribution to the g47g flux is not uniform over 
the entire field-of-view. The equivalent width of the emission lines 
in all three regions are given in TableQ] The central radio galaxy ex- 
hibits the emission lines with the largest equivalent width; Nv, Civ 
and Hell are prominent. Although galaxy A has a large Lya equiv- 
alent width it does not exhibit any metal emission lines, therefore 
emission lines would not contribute to the flux in the g475 pass- 
band. The diffuse continuum region shows strong Lya emission 
as expected since a large Lya halo encompasses the radio galaxy. 
Both Civ and Hell are detected with a combined equivalent width 
of 71 A, although there are variations in the strength of these lines 
along the slit length. As we are primarily interested in the nature 
of the intergalactic light (IGL) we adopt 71 A as the emission line 
contribution which means the emission lines are responsible for 17 
per cent of the total g475 flux. This contribution is subtracted from 
the g475 flux in the rest of this article. Since the contribution of the 
metal lines to the Spiderweb system halo is almost half that of the 
lMcCarthvl < fl993T) composite radio galaxy spectrum, we assume the 
contribution of metal lines within the Zg 14 filter to be negligible. The 
spectrum of the radio galaxy shown in the top panel of Fig.[2]shows 
numerous absorption lines redward of Lya. Because the interstellar 
absorption lines are generally stronger than stellar absorption lines, 
and none of the lines from the radio galaxy show P-Cygni like pro- 
files, we assume that these lines are interstellar and identify them as 
Sm/SiiiA1260, Oi/SinA1303, Cnll335, SiiiAl527 and A1A1671 
originating from gas at z = 2. 164. 



2.3 Extraction of galaxies 

Penteri cci et al. I dl998l) concluded that the UV bright clumps sur- 
rounding the radio galaxy were star forming galaxies. Defining 
which pixels belong mainly to a galaxy in this region is rather arbi- 
trary since the bright 'clumps' do not have regular shapes. We used 
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Figure 1. Left: Combined #4754-/814 image of the central 15"x 18.5" region of the Spiderweb complex. The radio galaxy, MRC 1 138-262, is marked with 
a white cross, and is surrounded by numerous satellite galaxies. The black box indicates the position of the MOS slit from which the spectra in Figure [2] 
are extracted. The label A indicates the galaxy A whose spectrum is shown in Figure [2] Right: §475 image of MRC 1138-262 and surrounding 15" x 18.5" 
region, smoothed with a gaussian of o~ = 0.05 arcsec. All regions identified as 'galaxies' by SEXTRACTOR have been removed and are coloured white. The 
intergalactic light is clearly visible between the radio and satellite galaxies. 




4500 
Wavelength A 



Figure 2. Rest-frame UV spectra of the radio galaxy, a nearby galaxy, and 
the region within the slit in which low level continuum and/or Lya emission 
was observed. The flux density of the diffuse continuum is per arcsec 2 . 
The wavelength range of £475 is approximately 4100— 5500A, therefore this 
band encompasses light from the Civ and Hell emission lines. 



Region 


size (arcsec 2 ) 


Lya 


NV 


Civ 


Hell 


Radio galaxy 


2.4 


147 


21 


57 


44 


Galaxy A 


2.4 


124 








Diffuse Continuum 


7.2 


202 




35 


41 



Table 1. Equivalent width of rest-frame UV emission lines in A from three 
separate regions: the radio galaxy, a nearby galaxy and the region in the 
slit in which low level continuum and/or Lya emission were observed. The 
sizes of the extracted regions are given in the second column. 



SEXTRACTOR dBertin & Arnouts|[l99r3) in the double image mode, 
where we use a detection image, constructed as the inverse vari- 
ance weighted average of the g 4 75 and /gi4 images, to determine 
the apertures, whilst photometry is carried out on the individual fil- 
ter images. A map of the total exposure time per pixel was used 
as the detection weight map. Regions with 5 or more connecting 
pixels which had a S/N greater than 3 times the local background 
RMS standard deviation were extracted. The right panel of Fig. Q] 
shows the 15"x 18.5" region surrounding the radio galaxy with all 
regions identified as 'galaxies' by SEXTRACTOR removed and col- 
ored white. AB Magnitudes of the galaxies were determined from 

hu = -2-5 log(counts) + 36.841 -0.074 (1) 
g475 = -2.5 \og(counts) + 36.856 -0.15, (2) 

dSirianni et al]|2005h . The last numbers on the right-hand side of 
both equations corrects for Galactic extinction. The slope of the UV 
continuum is defined as j8 given by = . Therefore a spectrum 
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that is flat in f v has = —2. We calculate the UV slope through 



-0.4 x (g 475 -7 814 ) 

lno ^475 



(3) 



where A475 and Asi4 are the effective bandpass wavelengths. Errors 
were calculated using the RMS images from Apsis and are quoted 
to 1(T unless noted otherwise. 

Since the galaxies were detected and extracted using SEX- 
TRACTOR isophotal apertures, only those pixels above a certain 
flux threshold were removed. Some of the flux in the wings of 
the galaxy profile is not removed by this process. As many of the 
galaxies have irregular profiles we cannot simply fit an ellipse to 
the galaxies and remove all flux within the ellipse. Instead we es- 
timate how much galaxy flux typically remains after the detected 
object is removed by S EXTRACTOR using three galaxies which 
have the most elliptical profiles compared to the other galaxies in 
the Spiderweb system. Azimuthally averaged radial profiles of the 
three galaxies were created, and compared to profiles of the same 
region after the galaxies were removed using our SEXTRACTOR- 
based algorithm. Fig. [3] shows an example of the radial profile of 
one galaxy (solid line) compared to the radial profile of the same 
region but with the galaxy masked out (dashed line). The dotted 
line indicates the continuum and background level. Fig. [3] clearly 
shows that some galaxy flux is not removed. The galaxy flux re- 
maining is estimated as the area under the dashed profile. On aver- 
age the remaining galaxy light totals 27 mag arcsec -2 in I$\4 and 
29 mag arc sec - 2 in £475. The galaxy's wings extend beyond the 
radius of the galaxy as defined by SEXTRACTOR across an area 
of approximately 0.6 arcsec 2 . Since the surface brightness of the 
galaxy's wings in £475 is 1.5 magnitudes below our definition of 
the detection limit of the intergalactic light (see section [2l4l . this 
remaining galaxy flux will not contribute greatly to the intergalac- 
tic light. However, the flux in the wings of the galaxies is brighter 
in 7gi4 compared to £475, therefore the §475-/314 colour can be af- 
fected: the regions surrounding the galaxies would appear to be 
redder due to the presence of the galaxy wings. No correction has 
been made for this effect, but the possibility of associated flux from 
the galaxies wings should be kept in mind when considering Figs. [4] 
and|5] 



2.4 Voronoi 2D binning 

To increase the signal-to-noise (S/N) the images were re-binned 
using the £475 image. The galaxies and the intergalactic regions 
were binned separately. The pixels of t he galaxies were binned us- 
ing the Voronoi 2D-binning method by Cappella ri & Copinl ( 2003) 
such that each bin had a S/N of 12. Many of the pixels at the cen- 
tre of the galaxies had a S/N much greater than 12, therefore they 
were not binned together with another pixel. The intergalactic re- 
gion had a lower S/N per pixel so the pixels were binned such that 
ea ch bin had a S/N of 5 us ing the Voronoi 2D-binning method 
by ICappellari & Copin (2003) with the Weighted Voro noi Tessel- 
lation modification proposed by Diehl & Statler (2006) instead of 
the Centroidal Voronoi Tessellation. We checked that there was no 
difference in any of the results if we binned the data using the /gi 4 
image. Using a much higher S/N (e.g. 10) resulted in very large bins 
which meant we lost a great deal of positional information. 

Since we wish to observe the intergalactic UV light we en- 
force that UV diffuse emission is detected in a bin only if the bin 
has a surface brightness brighter than 27 mag arcsec -2 in /g^and 
27.5 mag arcsec -2 in £475. This corresponds to a firm 5o" detection 
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Figure 3. /gj 4 radial surface brightness profile of the galaxy that lies at po- 
sition (2.3, 6.5) in Figs.|4]and|5] Solid line is the radial profile of the galaxy, 
whereas the dashed line is the radial profile of the same region with the 
galaxy removed by SEXTRACTOR. The dotted line indicates the continuum 
and background level. The area between the dashed and the dotted line is a 
measure of the galaxy flux that is not removed. 



of the IGL in the rest frame near-UV (NUV; 2215-3000A) and far- 
UV (FUV; 1300-1740A) bands above 5a confidence. In the rest of 
the article we refer to the IGL as the region in which the /814 flux 
is brighter than 27 mag arcsec -2 and the g475 flux is brighter than 
27.5 mag arcsec -2 , but not within a galaxy. 



3 RESULTS 

3.1 Distribution of the UV light 

SEXTRACTOR identifies 48 separate objects as possible galaxies, 
although only 25 are detected above the RMS noise. As we are 
interested in the light in between the galaxies we remove all 48 
detections of structure. The data in both bands were binned to a 
uniform S/N as described in section |2l4l and the surface brightness 
of the galaxies and the IGL are presented in Fig. [4] We define a de- 
tection of IGL as a bin which has 5(7 above the RMS background 
level in both g475 and 1^14. This corresponds to 27 mag arcsec -2 
in /814 and 27.5 mag arcsec -2 in £475. Because of this strong limit 
there may exist a fainter component that we do not have the depth 
in the images to unambiguously detect. In the 7gi4 there appears to 
be residual light around some of the galaxies therefore the stellar 
halos of these galaxies extend further than extracted by the SEX- 
TRACTOR algorithm (see section 12. 3t . Black areas in Fig. [4] are 
those in which the surface brightness is less than 27 mag arcsec -2 
in /814 or 27.5 mag arcsec -2 in g47g. The IGL surrounds the cen- 
tral radio galaxy in an ellipsoid shape, with best fit ratio 0.62 and 
angle 148.9°, centered on the core of the radio galaxy and extend- 
ing over 60 arcsec 2 (4243 kpc 2 ). Removing the galaxies within this 
region results in a total IGL area of 52 arcsec 2 (3168 kpc 2 ) The ori- 
entation of the intergalactic UV light is both similar to the Lya 
halo and within 30° of the radio axis. There is some extended light 
surrounding two other galaxies (one to the North, the other in the 
South, neither shown in Fig. |4), but generally the IGL is limited 
to the neighborhood of the radio galaxy and the group of central 
satellite galaxies. 

The total amount of continuum in the IGL is 21.7±0.1 mag in 
/§14 and 22.4±0.15 mag in g475. 39±7% of the total /g 14 -band light 
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Figure 4. 7gi4 surface brightness of the Spiderweb system. Colour bar scale is in mag arcsec -2 . Data has been binned in 2D so that each bin within a galaxy 
has a S/N greater than 10, and each bin outside a galaxy has a S/N of ~5 as described in section l2~4l Galaxies are outlined in black. Green contours are hya. 
The hya flux peaks on the galaxy that lies ~3 arcsec West of the radio galaxy. Image is rotated 5 degrees East from North, with North Up and East left. 

exists as the IGL, whilst 44±12% of the total g475-band light exists 
in the intergalactic region. 



3.2 Intergalactic UV colour 

Fig. [5] plots the UV slope /3 for all bins which have an I%u sur- 
face brightness greater than 27 mag arcsec -2 and a §475 surface 
brightness brighter than 27.5 mag arcsec -2 . It is difficult to distin- 
guish the galaxies from the IGL by their UV slope alone: the /gi4- 
8475 colour of the galaxy and the immediate area surrounding each 
galaxy are similar. For example the three red satellite galaxies to 
the East and Northeast of the central radio galaxy are surrounded 
by red IGL, whereas the blue galaxy to the West of the radio galaxy 
is surrounded by blue IGL. 

The average slope of the UV spectrum of the IGL is p g j = 
—0.84 ±0.5, but varies greatly from region-to-region. Fig.[5]plots 
the azimuthally averaged radial profile of the intergalactic UV light 
(all galaxies identified by SEXTRACTOR have been removed). It 
clearly shows that the colour of the IGL becomes bluer with in- 
creasing distance from the radio galaxy (marked by a white cross 
in Fig.Q). 



2 

Semi-major axis (arcsec) 



Figure 6. Azimuthally averaged radial profile of the slope of the UV spec- 
trum p. Errors are 3o~. The solid line has a slope of -0.37 and is the error- 
weighted best fitting line to the data points. The IGL becomes bluer toward 
larger radii. 
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Figure 5. Slope of UV spectrum (J3) of galaxies and intergalactic light of the Spiderweb system. Data has been binned in 2D as described in section Y2A\ 
Galaxies are outlined in black. The Lya flux peaks on the galaxy that lies ~3 arcsec West of the radio galaxy which has the bluest colour compared to the 
other satellite galaxies. Image is rotated 5 degrees East from North, with North Up and East left. 



Fig.[7]plots the UV continuum slope against the Ig\4 surface 
brightness. The crosses are the galaxies, and the diamond symbols 
are the averaged IGL grouped in bins of 0.25 mag arcsec -2 . The 
galaxies show a clear /3-surface brightness relation with a corre- 
lation co-efficient of —0.76. The jS-magnitude relation (not plot- 
ted) also shows a correlation but with a less significant co-efficient 
of -0.43. A similar relation is seen in a number of works (e.g 
Meur er et a j]|l999l :l Ouchi et "alll2004h . and possible causes may be 
a mass-metallicity, mass-extinction, or mass-age relation. This fig- 
ure shows that the galaxies are generally redder than the IGL. The 
average slope of the UV spectrum is jSg/ = —0.84 ±0.5, whereas 
60% of the galaxies in the Spiderweb system have a redder colour 
than this. Whilst the average colour of the total IGL is bluer than 
that of the average galaxy colour, the IGL directly surrounding a 
galaxy is often similar in colour to the galaxy itself. Therefore the 
IGL furthest away from a galaxy must have a bluer colour than the 
average galaxy colour. 



4 THE ORIGIN OF THE INTERGALACTIC LIGHT 

We have established the existence of low surface brightness UV 
light situated in between the galaxies of the Spiderweb system. 
Determining the origin of this IGL has parallels to the problem 
of the UV continuum excess of high redshift radio galaxies (see 
Tadhunter et al. 2002 and references therein). Numerous mecha- 
nisms can contributed to the UV excess including: scattered AGN 
light, young stellar populations, and nebular continuum emission. 
And the relati ve contribution of each mech anism varies between 
radio sources dSolorzano-Inarrea et al.ll2004r) . We now determine 
the origin of the IGL by investigating the possible emission mech- 
anisms. 



4.1 Faint cluster galaxies 

We examine the possibility that the intergalactic UV light may 
be caused by a fainter population of cluster galaxies that lie be- 
tween the detected satellite galaxies. Fig. [8] plots the rest-frame 
NUV (/814) luminosity function of the detected galaxies with the 
solid line showing the best fitting Schechter f unction 0(M) = 
O^log^lO ' 4 ^*-^]"- 1 exp[-10 a4 ( M *- M '] dSchechterl 19760 . 
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Figure 7. Crosses are galaxies, diamonds are IGL surface brightness 
grouped in bins of 0.25 mag arcsec~ 2 . The dashed line shows the detection 
limit imposed by the 27.5 mag depth limit in the §475 . The slope of the best 
fitting line to the galaxies is —0.9 ± 0.2 and the correlation co-efficient is 
—0.76. The best fitting line to the azimuthally averaged intergalactic light 
has a slope of —0.7 ±0.1. The slopes of both fits agree within lo~ of the 
error of the line fitting. Error bars on diamond points indicate lo" of the 
distribution in each surface brightness bin. 

The best fitting parameters were found to be M* = —24.3, <j>* = 
198.3 mag~'Mpc~ 3 and the slope a = —1.26. The slope of the 
faint end, a, i s found to be similar to that o f a sample of z — 2.2 
galaxies from ISawicki & Thompson! d20061) . The bright M* and 
large 0* are a direct consequence of the fact that the region of 
interest is associated with a massive galaxy at the centre of a 
forming cluster, and should not be taken as a measure of the lu- 
minosity function beyond this region. The luminosity function is 
complete to -19.5 mag and encompasses almost all of the galax- 
ies detected, therefore the integral of the Schechter function from 
— 19.5 to —24.5 results in the total 7gi4 flux from the galaxies in 
the Spiderweb complex. The intergalactic UV light in this band- 
pass accounts for 39% of the total flux above a surface brightness 
of 27.0 mag arcsec~ 2 . Therefore if the IGL resulted from a pop- 
ulation of fainter cluster galaxies, the missing galaxies would fill 
up the Schechter function up to —18.1 mag (indicated by the dot- 
ted line in Fig.[8]l. Therefore the IGL would comprise of approxi- 
mately 15 galaxies with absolute magnitudes between —19.5 mag 
and —18.1 mag. The IGL extends across a region 4243 kpc 2 (~ 
60arcsec 2 ). Therefore each of the 15 faint galaxies would have to 
extend across ~4arcsec 2 - larger than the massive radio galaxy, 
or be distributed in a very ordered arrangement around the radio 
galaxy. Such an ordered distribution is unlikely and any faint galax- 
ies are not likely to have such a large extent, therefore the IGL can- 
not be caused by a non-detected faint cluster galaxy population, but 
rather it must be due to an intrinsically diffuse source or very many 
unresolved radiation sources. 



4.2 Nebular continuum emission 

The distribution of the UV continuum and the Lya halo are sim- 
ilar which suggests that the continuum light may be produced by 
nebular continuum. Nebular continuum emission results from a 
combination of two-photon continuum emission, free-free emis- 
sion and free-bound recombination continuum, as well as a pseudo- 
continuum created by blended higher order Balmer lines just above 



Figure 8. The luminosity function of the galaxies. The solid line is a 
Schechter function fit to the data (see text for details). The dotted line in- 
dicates how much of the luminosity function is missing if the intergalactic 
light results from a population of fainter galaxies; approximately 15 galax- 
ies with absolute I magnitude (rest frame NUV) between —19.5 mag and 
— 18.1 mag would provide the necessary light. 

the 3646A jump. Theoretical case B nebular continuum emis- 
sion can be calculated u sing tabulated emission coefficients (e.g 
iBrown & Mathewslll970h and the measurement of the Hj3 emis- 
sion, therefore it can be well-constrained. We do not have a deep 
H/3 measurement of the Lya halo hence we estimate the average 
H/3 flux across the region from the measured Hell A 1640 flux in the 
diffuse continuum spectrum (Fig.[2]l. We use the STARLINK DIPSO 
command NEBCONT to estimate the contribution of the nebular 
continuum. Assuming a gas temperature of 15,000 K, density of 
100cm" 3 and 10% He abundance, the Hell flux of 8.46 x 10~ 18 
ergcm~ 2 s -1 arcsec~ 2 results in a nebular continuum in g 4 75 of 
~31.0magarcsec~ 2 or ~30.8 mag arcsec~ 2 in/814. Therefore less 
than 3% of the intergalactic UV light in both bands result from 
nebular continuum. 



4.3 Scattered light 

Previous studies have suggested that the UV excess in ra- 
dio g alaxies may arise from scattered quasar light (e.g Fabian 
Il989h and early polarimetric observations showed high degrees 
of linear polarization agreeing with the quasar scattering t heory 
dTadhunter et alll992tlcimatti et al.lll993l : IVernet et alfcOOll) . De- 
tailed studies of dimmer radio galaxies show that the contribution 
from scattered quasar light can be relatively minor dTadhunter et al.l 
2002). Here we examine the possibility that the diffuse IGL may re- 
sult from scattered light originating from either a hidden quasar or 
the stellar population of the radio galaxy as well as the numerous 
satellite galaxies. We first discuss whether it is energetically pos- 
sible for the IGL to result from either electron scattering or dust 
scattering of star light or a hidden quasar. 

The optical depth, t, of electron scattering is 

T = a T Rn e , (4) 

where o"x is the Thompson cross section (6.7 x 10~ 25 cm 2 ), R is the 
radial distance and n e is the electron density. The density of the 
gas in the Lya halo is esti mated to be ~100cm~ 3 but the volume 
filling factor is /,, ~ 10~ 5 dVillar-Martm et al 120031) . Therefore the 
average density over the halo volume is only 10~ 4 — 10~ 3 cm~ 3 . 
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However, any surrounding hot intergalactic medium (~ 1 x 10 7 K) 
may have a higher density (~ 0.15 cm~ 3 ) if we assume that it is 
in pressure equilibrium with the Lyof emitting material. Since we 
wish to determine whether electron scattering is possible we use 
the highest plausible value of n e in equation [4] i.e. that of the hot 
intergalactic medium. The light must be scattered up to 25 kpc from 
the central radio galaxy, therefore l~7.8x 10~ 3 . Assuming all the 
scattered light is observed as the IGL allows us to place a lower 
limit to the source intensity of 15.7 mag through 



IGL 



l-e- 



(5) 



The sum of all the UV light emitted from the galaxies (radio galaxy 
and the satellites) is only 1$ \a=1 1.2 mag, therefore the IGL cannot 
be produced by electron scattering of the stellar light from the ra- 
dio and satellite galaxies. Bright quasars at a redshift of 2 typically 
have an apparent magn itude of ~ 19, and the brightest is ~16mag 
(Schneider et al. 2007). Therefore it is energetically possible, al- 
though unlikely, that the IGL is produced by electron scattering of 
the light from a very bright hidden quasar. 

Dust has a far greater scattering efficiency than electrons, 
therefore if any dust exists in the halo, it will dominate the scat- 
tering process. Homogeneous optically thin dust scattering predicts 
a blueing of the emerging spectrum due to the wavelength depen- 
dance of the scattering efficiency of small dust grains. The IGL 
is observed to become progressively bluer with distance from the 
central radio galaxy (Fig. [7jl which ma y arise from multiple scatter- 
ing. Models from Cimatti et al. (1993) of dust-scattering of quasar 
UV light, where the scattered light is approximately proportional 
to the galactic light, ind icate a dust mass of ~ 10 8 M, ;i is required 
to scatter enough light. IVernet et al.l (2001) predict that only 7% 
of the incident light is scattered by 1O 8 M of dust in a spherical 
distribution of radius 40 kpc from the central source. Therefore a 
quasar of luminosity of 3 x 10 31 ergs -1 Hz~' (apparent magnitude 
of 19 at z = 2. 156) is required to illuminate the halo. It is plausible 
that MRC 1138-262 contains a quasar of apparent magnitude 19 
and submill i metre observations of both Stevens et al. (2003) and 
iGreve et alj 120061) show that MRC1 138-262 contains spatial ex- 
tended dust of total mass between 2 and 4.6 x 10 8 Mq. Therefore 
it is energetically possible that dust scattering of a quasar produces 
the IGL. Assuming a dust-to-gas ratio applicable to our Galaxy, this 
dust mass predicts a gas mass of 10 11 M Q , which is more than an 
order of magnitude greater than th e mass of ionized gas measured 
in this system: 2.3 - 6.5 x 10 9 M^ jNesvadba et alj|2006t) ). There- 
fore a cooler gas component must also exist in the halo. Evidence 
for cold gas is observed as absorption in the Lya emission lines of 
the radio galaxy, satellite galaxies and IGL shown in the spectra of 
Fig.H 

We now consider whether the model of dust scattering of 
quasar light is consistent with the morphology and colour of the 
IGL. A broad luminous Ha emission line with a full width at half 
maximum of 14, 900km s is a prominent featur e in the AGN 
spectrum of MRC 1 138-262 dNesvadba et alj|2006l) . however, the 
rest-frame UV spec trum (Fig. [2} reveals no nuclear broad emission 
lines. Additionally, Nesva dba et al.l d2006h measure a large extinc- 
tion (Ay = 8.7) from the luminous broad Balmer emission lines. 
This indicates that MRC 1138-262 harbours an obscured quasar. 
Therefore we expect any scattered light to be emitted from a bi- 
conical structure, revealing a butterfly morphology such as that 
observed in a number of low redshift Seyfe rt galaxies including 
the nearby NGC 5252 jTsvetanov et alJI 19961) . The smooth ellipti- 
cal morphology of the MRC 1138-262 IGL is not consistent with 



bi-conical scattering cones, therefore the morphology of the IGL 
argues against quasar scattered lig ht as a possible so urce for the 
light. In the clumpy torus model (Elit zuretal]|2004l) . where the 
torus is made of many small dense clouds, there is a finite probabil- 
ity that an unobscured line-of-sight to the AGN can occur through 
the torus. There have been a few recorded transits be tween type 1 
and type 2 AGN (e.g. Storchi-Bergm ann et al.lfl993"l) that may re- 
sult from such an opening in the torus. If such a situation occurred 
in MRC 1 138-262 the morphology of the scattered light would be 
elliptical dVillar-Martm et alj|2007l) as seen in the Spiderweb sys- 
tem. This scenario is unlikely as transitions from type 1 to type 2 
are rare, however this remains a possible explanation for the IGL. 

A further argument against the IGL resulting from quasar scat- 
tered light is provided by the uneven variation of the colour of the 
IGL. The IGL immediately surrounding a galaxy is very similar 
in colour to the galaxy (even though on average the IGL is bluer 
than the average galaxy colour). Such a clear correlation between 
the colour of the IGL and spatial position of the galaxies would 
require an unpleasant amount of fine-tuning if the IGL originates 
from quasar scattered light. Therefore the morphology and colour 
distribution of the IGL suggests that it does not originate from scat- 
tered quasar light. 

To summarise: the IGL cannot be produced by scattering (by 
either electrons or dust) of the stellar light of the radio and satellite 
galaxies as the scattering efficiency is too low. Scattering of quasar 
light is ruled out because MRC 1138-262 is believed to harbour a 
hidden quasar and the elliptical morphology of the IGL is not con- 
sistent with scattering from a hidden quasar. Furthermore the colour 
of the IGL shows a correlation with the colour of the satellite galax- 
ies, which is difficult to explain if the IGL resulted from scattered 
quasar light. We therefore conclude that the IGL does not result 
from scattered light from either a hidden quasar or the stellar light 
of the galaxy population. 

Scattering would result in high polarization, therefore imag- 
ing polarimetry measurements of the Spiderweb system would de- 
termine the importance and extent of any quasar scattered light. 



4.4 Synchrotron, inverse Compton, and synchrotron 
self-Compton emission 

4.4.1 Synchrotron emission 

The radio emission and spectral index can be used to determine 
whether the UV continuum emission may arise from the extrapo- 
lation of the synchrotron radio emission to U V wavelengths. The 
flux at 5 GHz is 1.6x 10~ 24 ergcm~ 2 s" 1 Hz" 1 dCarilli et alj 19971) 
and the integrated spectra l index is measured to be —1.8 from the 
4.2GHz to 8. 1GHz range jPentericci et all 19971) . This value is cor- 
roborated bv lCarilli et all 1 19971) . The core has a flatter spectral in- 
dex of —1.2 compared to the hotspots and extended emission. Ex- 
trapolating to 2600A, corresponding to the central wavelength of 
the 7gi4, results in a 14 flux of 35.0 mag. Additionally this mech- 
anism could not explain the IGL beyond the observed radio struc- 
ture. Therefore the UV continuum cannot be synchrotron emission. 



4.4.2 Inverse Compton 

Another possibility is the up-scattering of the cosmic microwave 
background photons (CMB) by the relativistic electrons of the ra- 
dio source. The energy density of the CMB increases at higher 



redshift by a factor of (1 



therefore inverse Compton emis- 



sion may be significantly more important at higher redshift. The 
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scattering electrons are presumed to belong to the same distribu- 
tion as those producing the synchrotron radio emission. Therefore 
the amount of UV flux Suv, at frequency Vuv, that can be pro- 
d uced through inverse Com pton can be estimated through eqn. 1 1 
of lHarris & Grindlavl fl979h : 

(5.05 x W 4 ) a C{a)G(a)(l+z) 3+a S R v% 



Suv = i 47 B l + « v g t/ =21.7mag, 

(6) 

where a is the absolute value of the spectral index (1.8), C(a) is a 
constant approximated to be 1.15x10 , G(a) is a slowly vary- 
ing function of a (G = 1.08 for a = 1.8), Sr is the radio flux 
at frequency v«, and B is the magnetic field strength. The mini- 
mum energy magnetic field strength B eq has been estimated though 
equipartition arguments to be ~ 200/iG JCarilliet all 1998b . There- 
fore it is energetically plausible that the UV extended emission can 
be produced by inverse Compton emission. It should be noted that 
the expected flux would be several orders of magnitude fainter if 
the spectral index is flatter. This up-scatting of the photons would 
require electrons with y ~ 1000. Electrons with y = 1000 would 
quickly lose energy through inverse Compton and synchrotron ra- 
diation in the environment of the IGL. It is difficult to devise a 
scenario where such a large volume is continuously populated by a 
diffuse medium of high-energy electrons. 

Although inverse Compton emission is energetically plausi- 
ble, both the morphology and colour of the IGL rule out inverse 
Compton as a plausible mechanism. The morphology of the IGL 
does not match the spatial distribution of the radio emission: the 
radio emission comes from the core, hotspots and jet-like regions, 
whereas the UV IGL has an ellipsoid shape centered on the radio 
galaxy and is brighter along lines connecting the satellite galax- 
ies. The radio morphology changes very little between 8.2GHz, 
4.7GHz, and 1.4GHz, and we assume that the structure of the radio 
source at low frequencies remains the same as that at high frequen- 
cies. Therefore we currently rule out inverse Compton emission 
based on morphological arguments, however, if a low frequency 
radio halo is discovered with the new low frequency radio observa- 
tories e.g. LOFAR, we must reconsider this conclusion. 

The most solid argument against inverse Compton emission 
relies on the slope of the UV spectrum of the IGL and its relation 
to the electron distribution index P defined as 



N{E)dE = CE~ p dE E x < E < E 2 



(7) 



where N(E) are the number of electrons with energy E between E\ 
and Ei- The spectral index (5) of the emission produced by single 
scattering of inverse Compton is related to the electron distribution 
index P by 

P- 1 



jB+2, 



(8) 



( Rvbi cki & Lightmanlll979h . and is identical to the relationship be- 
tween the spectral index (a) of the Synchrotron emission produced 
by the same population of electrons with distribution index P. S is 
trivially related to the observed slope of the UV spectrum, /3 . The 
observed range of /3 lies between —2 and 1, which corresponds to 
0<S<3orl</><7. Large variations in the colour of the IGL oc- 
cur across distances as small as a few kpc and such a large intrinsic 
range in P across these distances seem unlikely. Furthermore given 
the strong relation between a (or P) and the flux emitted in the 
UV (eqn. [6} any electron population which has P < 3 would con- 
tribute negligibly to the UV flux. Therefore regions of the IGL with 
j8 < — 1 cannot be produced by inverse Compton emission. Alter- 
natively we can imagine a situation where the IGL is produced by 



a population of electrons described by a single P and the range in 
produced from reddening by intrinsic dust. As a consequence of 
this we should expect to find a correlation between the UV flux and 
slope of the UV spectrum, such that the brighter regions would be 
bluer, but Fig.|7]shows the reverse correlation with the brighter re- 
gions having a redder colour. Therefore the colour of the IGL rules 
out inverse Compton emission as a possible origin of the IGL. 

4.4.3 Synchrotron self-Compton emission 

Synchrotron self-Compton emission occurs when the radio pho- 
tons are up-scattered by the same population of electrons that pro- 
duced the synchrotron radio emission. Although synchrotron self- 
Compton emission can also be ruled out on morphological grounds, 
we examine whether it is energetically possible. A condition for 
synchrotron self-compton to occur is that the ratio of the energy 
loss of an ultrarelativistic electron by inverse Compton to the syn- 
chrotron radiation (r\) must be greater than 1, i.e 

Urad Lrot 2/Iq 



'1 



> 1. 



(9) 



B 2 /2fl 4xr 2 c B 2 

U ra d is me radiation energy density, B is the magnetic field 
strength measured in t his hotspot through equipartion (10~ 3 Gauss, 
IPentericciet al. 1997), Hq is the permeability of free space, Lj ot 
is the total radio luminosity in the hotspot, and r is the size of 
the clump. We use the parameters of the UV continuum clump 
that lies ~3 arcsec West of the radio galaxy and matches the 
position of a radio hotspot. The radio luminosity Lj„, = 9 x 
10 44 ergs _1 jPentericci et alJT997b and r is estimated from the UV 
emission to be ~ l" 2 ~70kpc 2 . This results in rj < 1 x 10~ 4 , there- 
fore is not possible that the radio emission produced by the rela- 
tivistic photons are scattered to UV energies by the same popula- 
tion of electrons. 



4.5 Intergalactic star formation 

The final possibility is that the IGL results from UV emission from 
hot young stars. There are two possibilities for the origin of these 
stars: either they were stripped from the star forming satellite galax- 
ies, or the star formation occurred in-situ in the halo. The high 
rest-frame UV continuum suggests that we are observing OB stars, 
which have a main sequence lifetime of ~ lOMyrs and are thus 
good tracers of recent star formation. Stars recently stripped from 
galaxies would have a similar colour to that of the host galaxies, 
however the stripped stars would age passively whilst the galaxies 
continue to form stars, therefore the IGL should be redder than the 
galaxies in general if the IGL results from stripped stars. Since OB 
stars produce the vast majority of the UV emission and yet have 
a very brief lifetime, the colour of the UV light rapidly becomes 
redder once the star leaves the star forming region. Although the 
colours of the IGL immediately surrounding a galaxy is very simi- 
lar to the colour of the galaxy itself, generally the IGL is bluer than 
the galaxies as shown in Fig. [7] The virial velocity of a 10 M . 
galaxy, assuming a radius of lOkpc is 650kms~ 1 . An OB stripped 
star with a typical velocity of ~ 650 km s would travel a max- 
imum distance of 6.6 kpc before it turned off the main sequence, 
hence if the IGL were formed from stripped stars we would expect 
the regions furtherest from any galaxy to be reddest. Fig. [5] shows 
that the opposite is true: the regions furthest from a galaxy are bluer 
than regions close to a galaxy. This suggests that the stars that com- 
prise the IGL have not passively evolved, hence star formation oc- 
curs outside the galaxies in a > 60 kpc halo surrounding the central 
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radio galaxy, encompassing the Spiderweb complex. An alterna- 
tive possibility is that stars are stripped from the satellite galaxies 
together with molecular clouds which then form stars in the halo, 
however, this cannot be differentiated from star formation fueled by 
a gas reservoir in the halo. The similarity between the colors of the 
IGL just outside the galaxies and the galaxies themselves suggest 
some mixing of stellar populations or a common origin. 

There are a number of absorption features in the Lya lines 
shown in Fig.[2]which implies that cool dense gas coexists with the 
ionized gas in the halo surrounding the radio galaxy. These cool 
dense regions can be the fuel for the star formation. The rest-frame 
UV continuum flux at 1500A measured by the £475 band is a good 
tracer of star formation and can be converted to a star formation 
rate through 

SFR(M G yr- 1 ) = ^m_(ergs- 1 Hz- 1 ) (10) 

for a Salpeter initial mass function (Mada uetal.lll998t) . The total 
8475 band flux converts to a star formation rate of 130±13 M Q yr~' 
in an 4243 kpc 2 region with 57±8 M,:, yr~' of stars forming in the 
region between the galaxies (approximately 85 per cent of the area 
with a surface brightness > 27 mag arcsec~ 2 in /gi4, i.e. the total 
area of the galaxies and the IGL). Therefore almost half of the star 
formation traced by UV light occurs outside any galaxy-like struc- 
ture. 



4.6 /5-surface brightness relation for the galaxies and 
intergalactic light 

Having determined that the IGL originates from in-situ star for- 
mation we can now examine how p depends on metallicity, dust 
extinction and time since the starburst onset using the stel l ar pop - 
ulation synthesis models GALAXEV of iBruzual & Charlo 3 J2003h . 
We calculate /3 as we have defined it in equation[3] using the G and 
/ magnitudes outputted from GALAXEV. Fig.[9]displays the results 
of the GALAXEV modeling in which we determined the affects of 
variations in dust, metallicity and age. A star cluster formed at a 
radius of 30 kpc from the radio galaxy will be pulled by the grav- 
ity of the 10 12 Mq radio galaxy, reaching the centre of the halo at 
~ 5 x 10 7 yrs after it has formed (assuming it has no angular mo- 
mentum). This time is much greater than the lifetime of the OB 
stars which produce the UV radiation. Since the whole halo is UV 
bright the entire halo must be continuously forming stars. We there- 
fore use a model of a continuous star formation lasting for 1 Gyr. 
In all models we use the Chabrier initial mass function dChabrieJ 
l2003h . 

The left panel of Fig. [9] displays the relation between /3 and 
time since the onset of starburst. All plotted models have a metallic- 
ity of 0.4Zp,and three different extinctions are plotted: E(B — V) = 
0. 1 (solid line), E(B - V) = 0.25 (dashed line), and E(B - V) = 0.5 
(dotted line). The fainter dash-dot line marks the minimum time 
that a stellar cluster which formed on the outskirts of the halo would 
take to reach the centre. If the stellar cluster has any angular mo- 
mentum this time will be longer. The range in the UV slope that 
is observed may be produced by a constantly star forming stellar 
population if the stars at the centre are almost a Gyr old whilst the 
stars at the edge of the halo are only a Myr old, and the extinction 
is E(B — V) > 0.2. Therefore the gradient in /3 can be caused by 
an ageing of the stellar population if the star formation occurs on 
the outskirts of the halo, and the stars move inwards as they age. If 
the star formation rate has decreased over this period the range in p 



would be much larger than plotted in Fig. [9] since the older popu- 
lations would have a greater flux contribution from the older stellar 
population. Therefore the range in time may be less than 1 Gyr. The 
middle panel of Fig. [9] shows that the range of metallicity between 
0.2-1 Zq results in a change in j3 of ~ 0.4. Therefore metallicity 
variations in a young stellar population cannot produce the radial 
gradient in /3 for either the galaxies or the intergalactic light. The 
right-hand panel of Fig.|9]shows j3 verses extinction. A range in ex- 
tinction of AE(B — V) ~ 0.16 is enough to produce the range in jS 
observed for the IGL and just over double this amount will produce 
the entire range of /3 observed in the galaxies as well as the IGL. 

The IGL of the Spiderweb system is redder than an unob- 
scured starburst, the left-hand panel of Fig. [9] shows that a mini- 
mum E[B — V) ~ 0. 1 is required. Correcting the IGL for internal 
absorption assuming an E(B — V) ~ 0. 1 increases the IGL star for- 
mation rate to over 14OM0yr~'and more than 325M0,yr~' for 
the whole Spiderweb system. The range in observed j8 may be pro- 
duced by extinction or an aging stellar population, but the current 
data is unable to differentiate between these possibilities. 

4.7 Summary of results 

44%± 12% of the rest-frame FUV light within 8" (67 kpc) of the ra- 
dio galaxy exists as intergalactic light (IGL) that is situated between 
the radio galaxy and the satellite galaxies. In general the colour of 
the UV light of the galaxies is very similar to that of the surround- 
ing IGL, however, on the average colour of the IGL is bluer than 
the average colour of the galaxies, and there is a trend for the IGL 
to become bluer with radius from the radio galaxy. The source of 
this IGL is not faint cluster galaxies, but is either an intrinsically 
diffuse source or many unresolved sources that are spread out over 
a 3168 kpc 2 region. We have ruled out nebular continuum emission, 
synchrotron, inverse Compton, synchrotron self-Compton emission 
and scattering of either galactic stellar light or quasar light as pos- 
sible sources of the IGL. We conclude that the source of the IGL 
is most likely to be a distributed young stellar population situated 
between the radio galaxy and satellite galaxies, which has a UV- 
derived star formation rate of 57 M Q yr~' . 

Both the galaxies and the IGL show a /3-surface brightness re- 
lation which cannot be reproduced in stellar population models by a 
reasonable range in metallicity, but can be obtained by variations in 
either age or extinction, or a weighted combination of both. Either 
the halo is populated by a stellar population with a range of ages 
from a Myr at the outskirts of the halo to 1 Gyrs in the centre, and 
has an extinction of E(B — V) ~ 0.25. Or, irrespective of the age of 
the halo population, the extinction decreases from the centre to the 
outskirts by AE(B — V) ~ 0.16. Correcting for the minimum ex- 
tinction observed [E(B — V) ~ 0.1] increases the intergalactic star 
formation rate to over 140Me yr _1 . Most of the galaxies must ei- 
ther have a much higher extinction than this, or contain significantly 
older stellar populations. Both of these models can reproduce the 
radial profile of the UV color. 



5 DISCUSSION 

5.1 Relation between the halo of star formation and the Lya 
halo 

The spectrum of the Lya halo that encompasses the Spiderweb sys- 
tem contains strong Civ and Hell emission lines, which suggests 
that a harder ionizing source than stellar photoionization is at work 
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Figure 9. Range of UV slope j3 verses time since start of star burst (in units of 10 7 yrs), metallicity (in solar units) and extinction E(B — V). For first panel 
all lines have a constant metallicity of (1.4/ . . The solid line has an extinction E(B — V) = 0.1, the dashed line has E(B — V) = 0.25, and the dotted line has 
E(B — V) = 0.5. The vertical dot-dash line marks the minimum time that a stellar cluster which formed on the outskirts of the halo would take to reach the 
centre. The lines in the middle and right-hand panel have a constant extinction of E(B — V) = 0.2. The three different line-styles are cuts at different times 
since the start of the star formation. The solid line is lOMyrs, the dashed line is 200Myrs, and the dotted line is IGyr. The shaded region indicates the range 
of p observed in the majority of the IGL. 



here. However, the UV photons from the hot, young stars can con- 
tribute to the excitation of the halo and even dominate in some loca- 
tions. We can estimate the maximum Lya flux that can be produced 
by the total star formation rate of the galaxies and IGL by inverting 
the relation between Ha flux and the star formation rate 

Lffafergs" 1 ) = 1.27 x 10 41 SFR(M Q yr" 1 ) (11) 

dKennicutt et al.| [l994). Assuming case B recombination with no 
dust absorption ( Lya /Hq=1 1.4) and a total Lya luminosity of 
2.6 x 10 45 ergs -1 (Kurk 2003), this star formation rate can provide 
only 7% of the observed Lya emission. However the observations 
are unlikely to be complete and there may exist a fainter component 
of the IGL that has not been included in this calculation. Therefore 
7% may be considered a lower limit, although it is unlikely that this 
fraction can increase by more than a factor of 2. This leads to the 
conclusion that the Lya halo is powered by an additional mecha- 
nism to star formation. The spatial coincidence of a radio hotspot 
and the Lya peak has sparked speculation that shock heating may 
be involved. Other mechanisms for producing the Lya emission 
include AGN photoionization and gravitational cooling radiation. 

5.2 Low-redshift intracluster light and cD halos 

The low redshift intracluster light (ICL) is composed of old stars 
that are thought to have been tidally strippe d from merging gal ax- 
ies or the brightest cluster galaxy (BCG) jMihosetal.ll2005h or 
may have formed from Hi stripped from the galaxies during en- 
counters. The colour of the ICL is generally redder than the cluster 
galaxies and ther e is evidence for reddening of the ICL with radius 
from the BCG dKrick et al.ll2006h . The redder colour of the ICL 
agrees with the hypothesis that the stars were stripped from the 
galaxies and then passively evolved, whilst the galaxies continued 
to form stars. The diffuse UV light from MRC1 138-262 is gener- 
ally bluer than the satellite galaxies and the light becomes bluer 
with radius. This is inconsistent with the theory that the IGL stars 



were stripped from the satellite galaxies. It therefore appears that 
the these these phenomena have very different formation scenar- 
ios. The young stars in the halo of MRC 1138-262 cannot be the 
high-redshift analogues to the I CL observed at the cores of many 
low-redshift clusters, e.g. Virgo dMihos et al]|2005l) . However, the 
high-redshift halo stars may eventual become redder than the satel- 
lite galaxies if star formation ceases in the halo prior to the satellite 
galaxies. Therefore the high-redshift diffuse stars may contribute to 
the ICL at low redshift. 

Many BCGs show extended halos o f light beyond the de Vau- 
couleurs profile ( Matt hews et~al] 1 1964 ISchombertl 1 1988t) . These 
galaxies, such as M87, are known as cD galaxies. Deep obser- 
vations of M87 show that t he halo light gradually becomes bluer 
with radius dLiu et al.l20 05), similar to the IGL of MRC 1 138-262. 
However the observations of cD halos v ary widely in their conclu- 
sions. Both red and blue halos are found dSchomberJl98& Mackie 
1 19921), with red dGonzalez et al.ll2000T) and blue colour gradients 
dZibetti et alj|2005t) . In order for the halo of intergalactic stars to 
evolve into the halo of cD galaxies the halo must stop forming stars 
before a redshift of 1 and the population must a ge at least 5 Gyrs 
to reproduce the typical red colors of the halo. iLiu et alj d2005l) 
measure the envelope mass of M87 to be 3 x 10 11 M . Assuming 
the observed star formation rate of 80 Mq yr~ 1 is constant, it would 
take almost 4 Gyr to produce the necessary mass. Therefore there is 
just enough time for a cD envelope to form around MRC 1 138-262 
from the intergalactic star formation, then ages passively to have 
the same colour as other cD envelopes at low redshift. 

5.3 Origin of the gas reservoir: Cooling flows, cold accretion 
model, and outflows 

The high star formation rate of the halo requires that a large gas 
reservoir is present to provide fuel for the star formation. There 
are a number of theories as to how this gas reservoir may have 
arisen. As gas falls into the potential well of a dark matter halo, 
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it can be shock heated to the virial temperature, then radiatively 
cool, condense and form stars. This is th e conventiona l hot m ode 
of gas accretion. However, simulations of iFardal et aT] ( 1200 lh and 
iKeres et all d2005l) show that only half of the gas will be shock 
heated to the virial temperature of the potential well, whereas the 
rest will have a much cooler temperature (T< 10 5 K) as it accretes 
onto the galaxy. This mode of accretion i s referred to as th e cold 
accretion model (e.g lHaiman et alj|200d IFardal et alj|200ll) . Al- 
ternatively the large star formation rate of the radio and satellite 
galaxies can produce starburst winds, removing gas from the galax- 
ies and placing it in the halo. In this section we examine each theory 
in turn. 

Starburst mass outflow rates are rather uncertain, but var- 
ious observational techniques roughly agree that the mass out- 
flow rate is comparable to he star formation rate (see iHeckmanl 
120031 and references therein). Therefore a constant star formation 
rate of ~7OM0yr~' from the radio and satellite galaxies will 
result in an outflow rate of the same magnitude. The star for- 
mation efficiency of nearby molecular clouds is less than 10% 
( J0rge nsen et al .1120071) . therefore the current star formation cannot 
arise from a gas reservoir produced by a previous outflow unless 
the star formation rate was a factor of 10 larger at higher redshift. 

Cooling flows play an important role in the theory of galaxy 
formation. A cooling flow occurs when the shock heated infalling 
gas radiatively cools on a timescale longer than the dynamical time 
(Fabi anlll994l) . A prerequisite of a cooling flow is the presence of a 
hot thermalized medium. The high Faraday rotation measure of the 
polarized radio emission from MRC 1138-262 suggests that there 
is a hi gh column density wi th a magnetic field surrounding the 
galaxy jPentericci et alj 1997h . There have been X-ray o bservations 
of M RC 1138-262 with both ROSAT dCarilli etaljl998h and Chan- 
dra dCarilli et al.ll2002h satellites. Spatially extended X-ray emis- 
sion has confirmed an additional source other than the central AGN, 
however, there is no definite evidence for a thermalised medium: 
the extended X-ray emission from MRC 1 138-262 could have orig- 
inated from either inverse Compton emission or thermal emission 
from an ICM. 

A strong Lya emitting halo could be an observational sig- 
nature of the cold mode of accretion, and MRC 1138-262 is sur- 
rounded by a large Lya emitting halo. The star formation rate of 
the IGL and the galaxies is insufficient to power the Lya halo, 
therefore gravitational cooling radiation from the cold accretion 
may play a role. However, simulations have shown that the Lya 
radiation from gravitational cooling is less (or compar able in mas- 
sive s ystems) than that produced by star formation (Fardaletal. 
l200ll ), therefore cold accretion may be the origin of the halo, but 
we require an additional power source other than star formation 
and gravitational cooling radiation. Currently we cannot make any 
strong conclusions as to the origin of the gas reservoir, as all three 
proposed mechanisms cannot be ruled out. 

5.4 Metal enrichment of the intergalactic medium 

The data suggest that massive stars are forming up to 30 kpc away 
from the galaxy center. When these stars die they will release met- 
als into their surroundings, and because they are released so far 
away from the central galaxy the metals are much more likely to 
escape the gravitational potential and mix with the intergalactic 
medium. Numerical models show that high redshift massive galax- 
ies are i mportant contributors t o the enrichment of the intracluster 
medium jDe Lucia et alj|2004l) and inclusion of this radially dis- 
tant star formation may increase their influence. Further mixing 



may occur due to starburst winds. While the minimum average star 
formation rate density in the nebula is ~ 57MQVT /iz(30kpc) 2 
= 0.02 M^yr -1 kpc -2 , the star formation rate den sity may exceed 
the su perwind threshold of 0. lM Q yr _1 kpc~ 2 of iHeckman et al.l 
( 1990) when accounting for dust or in localized regions throughout 
the halo. 



5.5 Two modes of star formation in massive galaxies 

These observations show that star formation occurs in two modes 
in the high redshift universe: one in the usual Lyman break galaxy 
(LBG) clump-like mode on kpc scales, and the other in a diffuse 
mode over a large region surrounding massive growing galaxies. 
Almost 50% of the star formation traced by UV light occurs in this 
diffuse mode. There is a large amount of indirect evidence that sim- 
ilar star forming halos surround a number of other massive high- 
redshift galaxies: sub-mm measurements of a handful of z~ 2 radio 
galaxies show extended hot dust heated by star formation (Stevens 
et al. 2003; Greve et al. 2006); a number of z~ 2 radio galaxies 
have extended UV low surface brightness regions JPentericci et al.l 
1999) and slit spectroscopy of the Lya halos surrounding a large 
sample of high redshift radio galaxies show Lya/Hell r atios consis- 
tent with massive starbursts JVillar -Marti n et al 1 l2007l) . Lya halos 
similar in size and luminosity to the halos surrounding radio galax- 
ies have been found toward other high redshift sources a s well (e.g 
ISteidel et alj|2000l : iFrancis et alj|200ll; iDev et alj[2003) . It is be- 
lieved that these so-called Lya blobs are host to largely obscured, 
massive starb ursts with a possible contribution fro m superwinds 
and AGN (e.g lBower et alj2004l ; lcolbert et al]2006h . The new ob- 
servations of the Spiderweb system provide evidence that such ex- 
tended star formation can indeed occur in connection with massive 
forming galaxies and Lya halos. Long observations are required to 
detect the UV low surface brightness extended regions, therefore 
this phenomenon may occur frequently, but would easily be missed 
in previous observations. 

Clearly this mode of star formation does not continue uninhib- 
ited to low redshift as we do not observe much star formation in low 
redshift brightest cluster galaxies, especially outside the central few 
kpc. A number of mechanisms may halt star formation in the halo: 
the supply of gas may be exhausted, or a feedback mechanism from 
the star formation or the central AGN may blow away the reservoir 
of gas. At low redshift heating of the intracluster medium by the 
AGN is believed to be vitally important in pre venting a cooling flow 
dBinnev & Taboill995l : lchurazov et alj200ll) . In a related problem, 
some cosmological simulation s invoke AGN hea ting to halt star 
formation in massive galaxies dCroton et al]|2006l) . Although a ra- 
dio source is present in MRC 1 138-262, the observations of the IGL 
show that it is unable to prevent the gas from cooling and forming 
stars. The radio jets may not be able to couple to the surrounding 
medium as efficiently as it is supposed at low redshift, alternatively 
if the gas is accreting through a cold mode with a temperature of 
T< 10 K the radio phase of the AGN heating may not be powerful 
enough to prevent the gas from cooling, since the cooling function 
of gas at temperatures of 10 4 — 10 5 K is much stronger than at hotter 
temperatures. Additionally, the cold mode of gas accretion declines 
at low redshift, hence if the gas reservoir is fueled by cold gas ac- 
cretion, we only expect to observe such halos of star formation at 
high redshifts. 
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6 CONCLUSION 

The radio galaxy MRC 1138-262 is surrounded by a halo of 
IGL that extends across 60kpc. After examining nebular contin- 
uum emission, synchrotron, inverse Compton, synchrotron self- 
Compton emission, scattering, and stripping of stars as possible 
sources of the IGL, we conclude the most likely origin of the IGL is 
in-situ star formation. The minimum star formation rate (i.e. uncor- 
rected for dust) of the IGL is 57±8Mq yr~' , which is comparable 
to the total star formation rate derived from the UV luminosity in- 
tegrated over all the galaxies within 70 kpc of the radio galaxy (and 
including the radio galaxy itself). Applying a minimum dust cor- 
rection of E(B — V) ~ 0. 1 imply by the red colours of the the IGL 
increases the star formation rate of the IGL to 142M e yr _1 , and 
that of the whole Spiderweb system to more than325M Q yr- 1 .We 
estimate that the total observed star formation rate can produce ap- 
proximately ~7% of the Lya emission in the halo surrounding the 
galaxy. The radial colour gradient of the IGL indicates that there 
is a smooth range in extinction or stellar population age from the 
outer to the inner parts of the halo, where the inner parts are older 
or dustier than the outer region. 

While the presence of large quantities of ionized gas found 
around several remarkable species of high redshift galaxies has 
been attributed to a number of proposed mechanisms that include 
AGN feedback, infall, cooling flows, starburst superwinds and 
mergers, these observations of MRC 1138-262 show that any suc- 
cessful model should be able to accommodate the mode of extended 
star formation present in this paper. Our data suggest that the for- 
mation of the most massive galaxies is connected with that of their 
gaseous envelopes and star forming halos, part of which may pre- 
cede the intracluster light or cD envelopes, or perhaps may con- 
tribute to a satellite population. A significant amount of star forma- 
tion might be occurring in the form of extended low surface bright- 
ness features, beyond the typical UV detection limits, as well as in 
largely obscured extended halos as detected at infrared and sub-mm 
wavelengths. 
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